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Abstract 

The problem of a long, thin cylinder on a core, subjected to a finite 
band of axisymmetric radial load is considered. The core is analysed by 
elasticity theory. The cylinder is analysed by thin-shell theory. Appro- 
priate displacements and stresses are matched at the cylinder-core interface. 
After careful numerical evaluation of improper integrals, which appear in 
the solution, non-dimensional displacement and stress ratios for the core 
and cylinder are presented for different widths of uniform pressure band 
load. A comparison to existing solutions for the limiting case of a line 
load is made and shows the maximum radial stress to be significantly larger 
than previously reported. It is also found that an increase in the band 
width from the limiting case of a line load decreases the stresses and 
displacements in the assembly significantly even though the external resultant 
force in kept constant. 


Introduction 

In recent years the development of solid propellant rocket technology 
has been very rapid. Considerable attention has been given to the struc- 
tural safety and reliability of the motor case and its contents, which are 
a major component of such rockets. The present investigation is concerned 
with the effect of a uniform band of load intensity p and width 26 applied 
to a motor case assembly, as shown in Figure 1. Such a loading may arise 
when a clamp is placed over the assembly. The motor case is taken as a 
long, thin cylinder and its contents, consisting of solid propellant grain, 
is taken as a long solid cylinder or core. 

The core and cylinder are assumed to behave elastically. The core is 
examined by the use of equations from the theory of elasticity and the 
cylinder by the use of equations from classical thin shell theory. The 
elasticity equations for the core will be written in terms of displacements 
and solved with the aid of displacement functions. The shell theory equations 
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used to analyze the cylinder, are written in terms of its middle surface 
displacements in the radial and in the axial directions. It is then 
required that continuity of appropriate stresses and displacements at the 
interface between the core and cylinder be satisfied. In this way the 
cylinder unknowns may be related to the elasticity functions. 

The core displacement functions satisfy partial differential equations 
that will be solved by a separation of variables technique. The functions 
may be expressed in terms of coefficients of a Fourier integral representa- 
tion. When the load is also represented by a Fourier integral, the coeffi- 
cients may be evaluated from the shell theory equations which constitute 
the boundary conditions for the evaluation of the functions. 

To obtain desired stresses and displacements, improper integrals 
have to be evaluated. These are evaluated numerically after appropriate 
non-dimensionalization. A convergence test that established accuracies of 
one percent in the final results is incorporated in the computer program 
written for the numerical integration. 

In the limiting case as the width of band load becomes infinitely 
small, the assembly is loaded by an axisymmetric line load. The limiting 
case is of special interest and the results obtained are compared with 
those obtained by previous investigators. Comparison of the band load 
results obtained herein is also made to a plane strain solution for a 
finite cylinder. 

Field Equations of the Core 

The governing elasticity equations for the axisymmetric solid cylinder, 
or core, problem shown in Figure 1 are given below. For the cylindrical 
coordinate system (r,9,x) and W and U representing the radial and axial 
displacement components respectively, the strain-displacements equations 


3 . 


may be written 


3W „ W „ 3U „ 1,3W ^ 3U v ... 

£ r = 57 1 *B = r ; £ x = 37 5 £ rx = 2<37 + 57> (la_ld) 


where, £ , £ Q , £ and £ are the radial, circumferential, axial and 
* r 0 x rx ’ 

shear strain components respectively. The stress-strain relations are 


= Ae + 2Ge r ; ct q = Ae + 2G£g 


(2a, 2b) 


a = Ae + 2Ge ; T = 2Ge (2c, 2d! 

x x rx rx 

where, a , a n , 0 , and T are the radial, circumferential, axial and 
shear stress components respectively. The term e is three times the mean 
normal strain so that 

e = e + Eq + £ ( 3) 

r 10 x 

The terms G and A are the Lame material property constants. They are 
related to the modulus of elasticity, E c> and Poissons 's ratio, v, by the 
relations 

X - vE c 

A (l+v)(l-2v) 

G - £ c- 

2 ( 1+V ) (5: 

The stress components satisfy the equilibrium relations^ 1 ^ 

3 a r d T rx g 8- g r , c 

T - + — 5— =0 (6a. 

r 3 x r 


3 t rx 3 g x x rx 

+ + — = 0 


When the strain-displacement relations of Equations (1) and (3) are 
substituted into the stress-strain relations of Equations (21) one obtains 
the stress-displacement equations. 


G r 9W /o wW 3Uv, 

a = — 7- a.~+ (2-a)(- + 5—)] 

r a-1 3 r r 3 x 


*Numbers in superscripted brackets [ ] designate references in the Bibliography. 


KUWM f«Ui uitlrit •<>! >4 i WiH i 1 
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„ _ G r ,_ w 3W 9Uv W n 
9 “ a-1 [(2 a 9r + dx ) + °r^ 


„ G r / 0 w9W W, 9U n 

a = —rr L(2-a)(-r— + -) + a-r— ] 
x a-1 9r r 9x 


- r 9W 9U-, 
T rx = GC 97 + 97 ] 


(7b) 
(7c) 
( 7d ) 


where the quantity a is defined as 

a = 2(l-v) (8) 

Substitution of the stress-displacement relations of Equations (7) 
into the equilibrium relations of Equations (6) expresses the equilibrium 
equations in terms of the displacement components 


9 2 W 1 9W W a 1. 9 2 W ^ 1 9 2 U . 

Sr 2 ” r 9r r* a 9x^ a 3r9x 

2!» ♦ a . i)(£“ = o 

9x^ a^9r 2 r 9r a^9r9x r 9x' 


[ 2 ] 


(9a) 
(9 b) 


These equations were previously given by Marguerre , with a different 
notation, as his Equations (23'). 

Two functions (f x and f 2 ) of the coordinates r and x are now intro- 
duced in the following manner 


3f, 


9f 5 


W= r ^ + <2a ^> f 2 ] 


TT - 1 r 3fl ■ .. 3f> 

U * ' 2G C IT + r ix 3 


(10a) 

(10b) 


The introduction of displacement functions in this manner has been referred 

[3 4] 

to as a Two Function approach ’ . In these references it is shown that 
substitution of Equations (10) into (9) leads to the following relations that 
the functions f j and f 2 must satisfy 

9 2 f, 






(11a) 
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9 2 fa A 1 3 f> 1 


f 2 + Vx^ = 


(lib) 


The stress components may be written in terms of f j and f 2 by use of 
Equations (7). In view of Equations (10), the stress relations are 


a 

r 


a 2 fi 

r* 



+ a 


11 2 

3 r 


(12a) 


a = - llL 

9 r 3 r 


+ (1-a) 


IL 

3 r 


+ (1+a) 


1l 


(12b) 


a 

x 




(2-a) (^2. + — 

dr v 


(12c) 


T 


v>>. 


liL. 

3 r3x 


3 2 f, 


3 r3x 


+ (a-1) 



(12d) 


Equations (9) and (12) constitute the relationship of the displacer ent and 
stress components of the axisymmetric problem to the displacement .unctions 
fj and f 2 of the Two Function approach. 

Thin Wall Cylinder 

Classical thin-shell theory for a cylinder is often associated with 
the names of Timoshenko^^ and Donnell^ 6 ^. The derivation of the governing 
equations of this theory will be presented herein based on a simplification 
of the related problem as formulated by the theory of elasticity. This 
simplification is accomplished by the use of the following assumptions. 

In terms of geometry, stress, and displacement, the assumptions are made that 

(i) the wall thickness is considered negligible in comparison to the 
mean radius. 

(ii) the normal stress in the direction of small dimension is small as 
compared to the normal stresses in the other dimensions and may be 
neglected. 

(iii) plane sections initially normal to the undeformed middle surface 

remain unextended and plane and normal to the deformed middle * 

surface. 
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The elasticity Equations (1), (2), and (6) can be rewritten in terms 
of the shell theory notation for the z, y, x coordinate system shown in 
Figure 1, so that the strain-displacement relations are 


E yy ” R+z * E xx “ 3 x 


(13a-13c) 


zx 2 3 x 3 z 


(13d) 


The stress-strain relations are 


a = X e + 2G e 
zz s s zz 


0 = X e + 2G e 

yy s s yy 


(14a, 14b) 


0 = X e + 2G £ ; T = 2G e 

XX s S XX zx s zx 


(14c, 14d) 


while the equilibrium relations are 


3yzz . 9 Jy.v _ °yy - CT Z 
3 z 3 x R+z 


(15a) 


9 T zx 9 g xx T zx 
3 z 3 x R+z 


(15b) 


where, three times the mean normal strain e is related to the radial, 
circumferential, axial and shear strains respectively, by the relation 


e=e+e+£ (16) 

xx yy zz 

The quantities w and u are the radial and axial displacements respectively and 
z z 

R is the radius to the middle surface of the cylinder. The terms G s and X g 
are the material properties with the same definitions as given by Equations (4) 
and (5) where E is the modulus of elasticity of the cylinder and Poisson's 


ratio is given by y. 

Shell theory assumption (ii) indicates that O zz <* 0^ and O zz « 0 
so that 0 may be neglected. Rather than deal directly with the surviving 

ZZ 

stress components 0^, and t^, it is more convenient to deal with 
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resultants. They are defined as forces or moments per unit length of 


middle surface. N , N are the normal force resultants, Q is the shear 
x y x 

force resultant, while M^, are the bending moment resultants. They 

are related to the stress components by the definitions 

t t 


a dz 

yy 


N = 0(1 + £)dz ; N = 


= - T. (1 t £)dz 


o (1 + — )zdz ; M = - a zdz 

D * 1 ttir 


where, z is the radial coordinate measured radially outward from the mean 
surface of the cylinder and t is its wall thickness. The directions of 
these forces and bending moments are indicated in Figure 2. Shell theory 
assumption (i) states that t/R « 1 and since |z| £ y the quantities z/R 
may also be neglected in comparison to unity. Thus, Equations (17a), (18) 
and (19a) can be written simply as 


a dz 
_ xx 


Q x = - 


T dz 
. zx 


M = - I a zdz 


Shell theory assumption (iii) is generally called the Love-Kirchoff 


hypothesis for shells. From the assumption of inextensibility of the 
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cross-section, the radial displacement at a distance z from the middle 
surface is equal to the radial displacement at the middle surface. From 
the assumptions that plane sections remain plane and that sections initially 
normal remain normal, the axial displacement at a distance z from the middle 
surface can be related to the displacements at the middle surface. Denoting 
the radial and axial displacements at the middle surface by w ard u respectively, 
one concludes that 

d w 

w = w ; u=u-Zt— (23a, 23b) 


3w 

u z = u ' 2 37 


Sutstitution of these equations into Equations (13) and neglecting z com- 
pared to R in Equation (13b) by shell assumption (i) the strain components 
become 

e__ = 0 ; e__ = £ (24a, 24b 

zz 


w 

"yy ’ R 


3u z3 z w < y 

xx " 3x 3x^ * zx 


(24a, 24b) 
(24c,24d) 


Since it has been assumed that a « a , and a « a so that 

zz xx zz yy 

a vanishes, the constitutive elasticity Equations 04b) and (14c) together 

ZZ 

with Equation (14a) lead to 


a = -r-— r [e „ + ye ] (25a) 

xx 1 ~jj 2 xx yy 

ct = : ■ E - * [e + ye ] (25b) 

yy 1 -m 2 yy H xx J 

The stress-strain equations which involve z quantities (Equations (14a), 
(14d)) can be seen to be violated. This is an expected consequence of the 
additional assumptions that have been made for shell theory for quantities 
in the z direction. Substitution of the strain-middle surface displacement 
relations of Equations (24b), (24c) into the stress-strain relations of 
Equations (25) givas the stress-middle surface displacement equations 

E r 9u z3 2 w . uw, 


a = E r 2H. - + H] 

a xx 1-y 5 l 3x Sx 5 R J 


(26a) 
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E_ rW y3u pz3 2 w 

yy ' 1-ya L R 3x ~ 3x2 J 


(26b) 


Substitution of the foregoing into Equations (17b), (19b), (20) and (22) 
which define the stress result: s,and performance of the indicated inte- 

grations, shows that the stress resultants are related to the middle surface 


displacements as 


N 


x 


Et 

1-M 2 


r 3u pw-, . M 
C 3^ + "R ] • N y 


Et r p9u w-, 

l-y2 L 3x R J 


M 

x 



M 

y 


D 3 2 w 


(27a, 27b) 
(27c, 27d) 


where, the flexuial rigidity of the shell D is defined 


D = 


Et 3 

12(1 - pO 


(28) 


Equilibrium requirements in the x and z directions as well as moment 
equilibrium about an axis in the y direction, as shown in Figure 2, demands 



+ P. 


0 


(29a) 




0 


(29b) 



(29c) 


where, p and p are tbe radial and axial loads per unit middle surface area, 
z z 

These are shown in Figure 2. Substitution of Equation (29c) into Equation (29b) 
gives 



(30) 


Substitution of Equations (27a-27c) into Equations (29a) and (30) gives 


R 2 3 2 u 

5x2 


+ HS&+ (iV) £ 


3x 


Et 


p = 0 

r v 


Ola) 


yR3u ^ R 2 t 2 3 l *w 
^ — + w + -rr- -5 — ir 
3x 12 3x* 



z 


0 


(31b) 
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When no axial surface load is present (i.e. p x = 0) these equations reduce 
to the single familiar equations for axisymmetric cylindrical shell problems 
and beam on elastic foundation problems. This is shown and used in Appendix 1. 
Boundary Conditions 

Equations (10) and (12) for the core and Equations (26) and (31) for the 
cylinder are the relations that must be satisfied to obtain solutions for the 
stresses and displacements. The core unknowns have been expressed in terms 
of the displacement functions fj and f 2 . The boundary conditions provide 
the means of relating the cylinder unknowns to these functions and necessary 
relations for the evaluation of the functions themselves. 

The following boundary conditions for the cylinder-core problem must 
be satisfied. The stresses and displacements must all remain finite as 
x -*■ ± ”. Continuity of radial and shear stress as well as continuity of 
radial and axial displacements must exist at the cylinder- core interface 
for all values of axial coordinate. Let the core radial and shear stresses 
at the outer radius of the core (r = a) be designated a and T and let the 

3 3.X 

radial and axial displacements at the same surface be designated U and W . 

3 3 

The stresses at the cylinder-core interface are assumed to act on the cylin- 
der middle surface. For an axisymmetric externally applied radial load or. 
the cylinder p*(x), chosen outward positive in accord with the coordinate 
system, the cylinder unit loads are 

p = p*(x) - a ; p = -T (32a, 32b) 

* a *x ax 

From Equations (23) the cylinder displacements at the interface are given by 


w 

z 


z 




t 3w 

u + 2 


(33a, 33b) 


Since the left hand side of Equations (33) must equal W and U respectively, 

3 3 
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by the displacement continuity conditions above, the middle surface dis- 
placements of the cylinder may be written as 

“ = “a > u = (34a,34b) 

When the condition imposed by Equation (34b) is used the results obtained 
are referred to as an interface matching solution. 

It is in the boundary condition expressed by Equation (34a) that pre- 

[7 gT 

vious investigators of the line load solution 1- * assumed 

u = U (35) 

a 

When this boundary condition is used the results obtained are referred to 
as a middle surface matching solution. The axisymmetric line loading of a 
cylinder with no core helps to shed some light on the difference in the 
matching. In Figure 3 the axial variations of the axial displacement para- 
meter 2Gu z /q x 10 3 are plotted for the cylinder inner and middle surfaces. 
The quantity q is the resultant force per unit circumferential length. The 
results in this figure were obtained from shell theory expressions derived 
in Appendix 1. The effect of the difference in the use of Equation (34b) or 
(35) in the final results of the cylinder-core combination will be evidenced 
once numerical results are obtained. 

In addition to the boundary conditions for the stresses and displace- 
ments it is necessary to relate the radii of the outer radius of the core, 
a, and the radius to the middle surface, R. Since 

a = R (1 - ~) (36) 

and it is assumed that t/R « 1 for thin-shell theory it is taken that 

a = R (37) 

The cylinder displacements can be expressed in terms of the elasticity 
stress functions fj and f 2 through the substitution of Equations (10) and 
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(34) into Equations (33). No additional information for the redial displace- 
ment would result since this can be evaluated as the core radial displacement 
at the outer core radius. When the following definitions are made 


f ia = f > 


f = f 
2a 2 


(38a, 38b) 


the axial displacement is given as 


u = - ~ [f + af - (z + -i-Hf + af - ( 2 a-l) f )] 

z 2G ia,x 2 a, x 2 ia,rx 2 a, rx 2 a, x 


The cylinder stresses can similarly be expressed in terms of f ja 
and f defined in Equations (38). By the substitution of Equations (10), 
(34), (37) and (39) into Equations (26), it is found that 

a = G {-[f +af ] + (z + |)Cf +af -(2c-l)f ] 

xx m la, xx aa,xx 2 ia,rxx 2 a,rxx 2 a, xx 

-U Ef, + af -(2a-l)f ]} (40a) 

— ia,r 2 a, r 2 a 

u 

a = G { - — [f + af -(2a-l)f ] -y[f + af ] 

yy m a ia,r 2 a, r 2 a ia,xx 2 a, xx 


ni( z + ^-)[f 


2 ia,rxx 2 a,rxx 


-( 2 a-l)f„ ]} 


2 a, xx' 


(40b) 


where the quantity G m is defined by 


Equations (39) and (40) thus express the cylinder axial displacement 
and stresses in terms of the elasticity functions fj and f 2 . 

Unknown Coefficients and Their Evaluation for a Band Load 

The solution for the core unknowns can be achieved with the aid of the 


functions f Jt f 2 which must satisfy Equations (11). These functions can be 

[3 4 ] 

solved for by a separation of variables technique * . If even functions of 
the axial coordinate are considered the expressions for fj and f 2 can be written 
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A I (n — ) cos (n *) dn 
o a a 


f 2 s 


r b ii 
1 0 


(n — ) cos (n —) dn 

S 3 


(42a) 


(42b) 


where A and B are functions of n and I and I. are modified Bessel Functions 

o 1 

of the first kind or order one. The quantities A, B depend upon the continuity 
relations between the core and cylinder. Before the evaluation of these 
quantities is considered, it is beneficial to express all the displacements 
and stresses of interest in terms of A and B. 

For the core, substitution of Equations (42) into the displacement an 
stress relations of Equations (10) and (12) and noting the derivatives of 
the Bessel functions 


3[I (n -)] 


o a _ n 


= ii i ( n L.) 

3r a 1 a 


(43a) 


3[Ii (n ~)] 


■ ■ *-■- = — I (n -)-■=■ I,(n -) (43b) 

3r a o a r 1 a 

it i; found that the radial and axial displacements of the core are given 

by the expressions 


2GU 




[AI (hT) + B— I (n^*)] sin (i£) dn 
o “ a i a a 


(44a) 


= f {-AnI l 
* n 1 


2GW = I {-AnI (n^-) + B [2a I (r£) -i£ I (n^)]} cos (r£) dn (44b) 
I 1 a iaaoa a 


The radial, circumferential, axial and shear stresses within the core are 
given by the relations 


= f(A [ - 4 I (r£) + £ Ij <«£)] 
I ci o 3 r * a 

• n 


-B [(^-*r + “-)Ii (n^)-(l+o) ~ J (n^)]} cos (nr - ) dn 


a o a 


(45a) 
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{-A— I (i£) + B[(l-a) — I (n—) + -7 I,(nJ)]} cos (i£) dn (45b) 
r 1 a aoa r 1 a a 

= f{ A— I <i£) + B[(2-0) £ I (n^) + I (nJ)]}co: 
x aoa aoa a 2 1 a 

J 0 

(A— I.(r£) -B [c£ I.(r£) -n 2 I (nj)]} sin (nj) dn 
a 1 a a 1 a a oa a 


A 

a 


(45d) 


It can be noted that similar expressions for the displacements and the 

[7] 

shear stress were given by Yao for the line load solution. For the line 
load solution particular values of the coefficients A and B apply. A dis- 
crepancy in the radial stress expression he used exists. This was carried 
through by Yao in his subsequent equations and calculations. 

Since the quantities U , W , o and T are necessary for the boundary 
conditions expressed by Equations (32) and (34) they may be written simply as 


■r* 

* n 


2GU = I n[AI (n) + BI (n)] sin (n— ) dn 
a j o 1 a 


(46a) 


2GW = 
a 


{-Anl,(n) + B[2al (n) - nl (n)]} cos (r£) dn 
1 *0 a 


(46b) 


a =i 

a a 


{A [nl,(n) - nl (n)] 
1 o 


-B [ (n 2 +2a) I (n) - (l+a')nl (n)]} cos (i£) dn (46c) 

1 o a 


T 

ax 



(An 2 I (n) ~B[otnI (n) - n 2 I (n)]} sin (r£) dn 
1 1 o a 


(46d) 


The cylinder axial displacement and stress expressions in terms of 
A and B can be arrived at by the substitution of Equations (42) into 
Equations (39) and (40). This gives 
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U z = 2G 


{A [nl (n) -nl,(n)3 

O * 


+B [nl.(n) +n(2ot - I (n))]} sin (r£) dn (47a) 

* no a 


a = G , 
xx m j 


(a [ i (n) - r (y+n) I (n)3 


a o 


-B [ — (y+n) I (n) - — (n 2 +2a(y+n)) Ii(n)3) cos (n-) dn (47b) 


a = G 

yy m 


o 

{A [ Hsi I (n) - f (l+yn)I. (n)3 
a o a * 

-B [ — (l+yn)I (n) - — (yn 2 +2a(l+yn))I,(n)3) cos (r£) dn 


(47c) 


where the quantity n is given by 


n 


n 2 (2^+ 1) 


2(f) 


(48) 


From Equations (44), (45) and (47) and a mathematical theorem pre- 

[93 

sented in Sneddon the displacements and stresses can be shown to satisfy 
the boundary condition requirement that they remain finite as the axial 
coordinate goes to infinity. 

The radial load, which up to this point could be any general axisymmetric 
distribution, can be represented by a Fourier integral. In correspondence 
with the forms of the expressions for fj and f 2 given by Equations (42) the 
band load representation is taken as 


p*(x) 


JL 

ira 


#00 

Ceos (n— ) dn 

* n 3 


(49) 


where, C is a function of n given by 



(50) 


This expression for the Fourier integral representation of an even function, 
to which the load is restricted, may be found in any mathematics text in 
which Fourier integrals are discussed (e.g. Kreyszig^ 10 ^ ) . For the particular 
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case of a band load of uniform 

C 


where q is the resultant force 


Figure 1 and 26 is the total axial width of the load with 


q = 26p 

(52) 

Substitution of Equations (32), (34), (37), (46), (49) 

and (51) into 

Equations (31) leads to two integral relations containing 

A and B, namely 

[ " 1g ^n A + ^i 2 B ^ dn = 0 

* 0 

(53a) 

j [X^A + A 22 B - ?] cos(n-) dn = 0 

j qZo a 

(53b) 

where the parameters 


A u = -n 2 I Q (n) + (g + ~^)n Ij (n) 

(54a) 

X 12 = (g + ~~) n I Q (n) - Cn 2 (l+ —•) + a(y+g)J I^n) 

(54b) 

X 2J = n 2 g I Q (n) - [1+ + kn 4 + g-y]n I^n) 

2 

X 22 = -[1 + + kn 4 -(l+a)(y-g)]n I Q (n) 

(54c) 


+ [(n 2 +2a)g + 2a (1-y+kn 4 ) + - n ° * ]I.(n) (54d) 

a * 


and the material and geometric parameters C* k and g are given by 



(55a) 

k = ii?r 

(55b) 

g ■ y -2(i-y 2 ) | 

(55c) 


pressure p, C reduces to 

sin(n^) 

= -q - V- 

(n I ) 


(51) 


per unit circumferential length as shown in 


If it is taken that the axial displacement matching condition is the 
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middle surface matching given by Equation (35) rather than Equation (34b) 
some of the terms of Equations (54) would not apply, but instead one would 
find 

A u = -n l 2 * I Q (n) + gn Ij(n) (56a) 

X 12 = gnI Q (n) -[n 2 +a(y+g)] I x (n) (56b) 

A 21 = n 2 g I Q (n) -[1 + kn**+ g-y]n Ij(n) (56c) 

A 22 = "Cl + kn 4 + (l+a)(p-g)]n I Q (n) 

+ C(n 2 +2a)g + 2a(l-y+kn 4 )] Ij(n) (56d) 

To satisfy Equations (53) it is necessary that 


A U A + A 12 B = 0 (57a) 

A A + A,,B = Z, (57b) 

21 zz 

The direct evaluation of the coefficients A and B which are functions of 
n is not necessary. Substitution of the algebraic solution of the simple 
linear relations of Equations (57) shows that 


- _ 1 £ 


where the denominator, is given by 


d a 


(58a, 58b) 


°A = 


ll 


'12 


= ^ 1 1^22 " ^ 12^21 


(59) 


l 21 *22 

Appropriate values of A and B can be substituted into the desired displace- 
ment and stress expressions. When the values of A and B given by Equations 

(59) are substituted into Equations (44), (45), and (47) the final forms of 
the expression to be evaluated numerically result . These are easily obtained 
and are all included here in nondimens icnal form for completeness. 

The expressions for core displacements are given in nondimens ional form 


as 


2GU 

q 


1H 


sin(n^) 

a 


dn 


(60a) 
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* fw 4 cos (n£) dn 
q J.nD» a 


* *o “ "X 
where the terms under the integral are 

U n 3 nC “ X ia I 0 (,, J ) + X n 7 I » (n 7 ):i 


W s X n I, (»r) + X . . [ 2a I (n— ) - r>— I (n—)3 
n 1 2 *i i a a o a 


and the parameter C is defined 


F = i 


The corresponding nondimensional core stress relations are 


afli 

q 


* j (a ) 5- cos(i£) dn 
Jo r n D X a 


q 


* r<v„ 


cos (nr) dn 


aa» 

q 


■ l>.’. 


4- cos (nr) dn 


* f°(T ) sin(n^) dn 
q J n rx n D, a 


'o J X 

where the terms under the integral are 

(a ) = X. _[n 2 I (i£) - I.(n^)] 

r n 12 o a r/a * a 


•X.-C^-t—) I.(nJ)-(l+a)nI (nj)! 
a r » a oa 


<V„ ■ Krthty * x » C( 1 -° ,nI o <n i ) + r 


(ff. 


) = -X,,n 2 * (»7> ♦ X-.C^-^nl (x£) + n 2 J I x (i£>] 

t n 12 o a ii oa a 1 a 


<*„>„ * - x i,n ! I,<>T) * X u CcmI i (A - n*f 


(60b) 

(61a) 

(61b) 

(62) 

(63a) 

(63b) 

(63c) 

(63d) 

(64a) 

(64b) 

(64c) 


(64d) 
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The computed core values to be presented herein are to be evaluated 
at the interface, r = a, only. For this case, Equations (61) for the 
displacements and Equations (64) for the stresses reduce to 


U 

n 


n[-X 12 I Q (n) t X ll l J (n)3 


r=a 


= A 


r=a 


(a) 

r n 


r=a 


,,nl (n) + A [2al (n) -nl (n)] 

1 2 i 11 1 O 

= A [n 2 I (n) - nl (n)] 

1*0 i 

-A [(n 2 +2ot)I (n)-(l+a)nl (n)] 
11 l o 


(65a) 


(65b) 


(65c) 


(0 e ) n 


(o ) 

x n 


= A.,nl.(n) + A.,[(l-a)nl (n) t 2al (n)] 

12* 11 O 1 


r=a 


r=a 


-A,,n 2 I (n) + A,.[(2- 0l nl (n) + n 2 I.(n)] 

12 O 14 O 1 


(65d) 


(65e) 


(T ) 
rx n 


= -A n 2 I (n) - A [anl.(n) - n 2 I (n)] 

i*i ill .i 


(65f ) 


r=a 


The nondimensional cylinder axial displacement and stresses are 
given by 


2Gu z _ 

l (u } D 

sin(n— ) dn 

(66a) 

q 

Jo 2nD A 

a 


an 

- SX = 

G f°°(a ) 

cos(n^) dn 

(66b) 

q 

ml xx n 

1 0 

D x 


a o 

- 

G T(a ) 

I*— cosfrr— ) dn 

(66c) 

q 

m J yy n 

D. a 



where the terms under the integral are 

(u ) = -A Cnl (n) - nl (n)] + A-.Cnl^n) + - I (n)J (67a) 

zn 120 l hi n o 


( a xx^n = “^i2^ n2l o^ n ^ " n ^+n)I l (n)] 


-A u Cn(p+n) I Q (n) - (n 2 +2o(y+n))Ij (n)] (67b) 
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(a ) = (n) - n(l-Hir| )I, (n)] 

yy r * 2 o 1 

•X u Cn(l+yn)I o (n) -(yn 2 tfc(l+vn ) )Ij (n)] (67c) 

The computed cylinder values to be presented herein are the stresses which 

2z 

are evaluated at the cylinder outer surface, — = 1 , only. For this, the 
quantity n from Equation (48) reduces to 


n 


2z 

t 



( 68 ) 


Computational Difficulties 

The formal solution to the band of load on a cylinder and core assembly 
as derived in the previous section is sufficiently complicated to cause 
computational difficulties. The methods used to overcome these difficulties 
will be discussed in the present section. 

It can be seen from Equations '60), (63) and (66) that a typical inte- 
gral requiring solution may be expressed as the nondimensional quantity S 
defined as 


■r 

J n 


s.in (n^-) sin (i£) 

s{(^“) t — — ){ or a } dn 

o (n^) cos (i£) 

3 cl 


(69) 


where, s is the quotient of two expressions that involve powers and modified 
Bessel functions of n. It may be written as 

s_ 


s = 


(70a) 


where s r takes on values given by Equations (61), (64), and (67). General 

computational difficulties encountered in elasticity solutions for cylinders 

[3] 

were discussed in a previous report by the authors*" . In that report, it 
was indicated how the modified Bessel functions can be evaluated and the 
general method to be used herein for the accurate evaluation of the improper 
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integrals was given. At this point the details of the numerical evaluation 

of the particular typical integral of the cylinder-core problem in Equation 

(59) are given. The difficulties that arise are caused by the function s 

and by the presence of the trigonometric functions for both band width and 

for axial location. They will be discussed in this order. 

For both 'small' and 'large' values of Bessel function arguments, 

values of the numerator, s and denominator, D, of s too small or too 

n A 

large to be handled by a digital computer will occur. To elminate this 
difficulty both numerator and denominator in s are divided by n 3 I Q 2 (n). 
Equation (70a) is rewritten to read 



(70b) 


where, -he nrmerator and denominator are given by 


S n ~ n 4 I 2 (n) 
o 


n J I 2 (n) = * 11*22 ' ^ 12^21 


The A's are taken as 


r u , , 2 t 

X U = - 1 t(g + n 2? )g 


12 


12 - nTTn) = <8 = ™<u + g)]i 


o 

- = * 21 
2 i " n 2 I (n) 
o 

A 22 


2 

= g -Cl + C2y-g) + kn 1 * + g-y]g 


n 2 t 


l 22 * nTTnT = - [1 + (2M “ 3) + kn4 ’ (1+0t)( ^ )] 


+ [(n 2 +2a) g + 2a(l-y+kn 1 *) + ~~ (3u-g)]g 


where, g is the ratio of modified Bessel functions given by 

g 


Ij (n) 


nl (n) 
o 


(71) 

(72) 

(73a) 

(73b) 

(73c) 


(73d) 


(74) 


Again some terms would have been lost had displacements been matched 
with the middle surface so that the A's would be given as 
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*11 

n " n*ljn) = ' 1,g ‘ 


.2 = Ki^T = B - Cn ' + “ (u + g)] * 

X. 


21 


2 i “ n z I Q (n) 


= g -[1 + kn 4 + g - y]g 


{ 22 ~ nl (n) ^ + kn (1 + a) (U g)] 

o 

+[(n 4 + 2a) g + 2a (l-p+kn 4 )]g 


(75a) 

(75b) 

(75c) 

(75d) 


The values of the displacements necessary for computational purposes 
for small and large arguments from Equations (61), (71) and (73) are 

(76a) 


“n = * \. <n t )2 * ] 

"n = - io ]} 

where g Q and g x are the ratios of modified Bessel functions given by 


(76b) 




I (r£) 

0 a 

1 (n) 
o 


g i = 




(n— ) I (n) 
a o 


(77a, 77b) 


Similarly for the core stresses, from Equations (64), (71) and (73), the 
values for small and large arguments are 


(8 r ) n * A 12 [g 0 - i,3 C(n|) 2 ♦ 20] g, -U + o) g o ) (78a) 

< 3 e>„ = + - “> * 0 + 2o «. ] (78b) 

<Vn = - 3 ‘ 1 26 o + *„«*-«»„ + (r f )2 S. ] <78C) 

(f rxV "I C ‘ X i2*l ‘ (76d) 

At the outer radius (r = a) g Q = 1 and g t = g. The quantity g, given 

by Equation (74), is evaluated for small values of n (less than .001) by the 

division of Equation (42b) of by Equation (42a) of For large values 

[3] 

of n (greater than 50) g is evaluated by the division of Equation (45b) of 


[3] 


For ^ < 1 and n < .001 the quantity g o 


by n times Equation (45a) of 
may be evaluated by the divis 
same Equation for H = n. For the same condition g i may be evaluated by the 


[3] — r 

may be evaluated by the division of Equation (42a) of for r) = n^- by the 
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division of Equation (45b) of^ 3 ^ for n=rr- by Equation (45a) of^ 3 ^ for n=n. 
For r/a < 1 and n > 50 two methods of evaluation of each of g Q and gj 
must be considered. These methods depend upon the magnitude of n r/a. 

For n^- £ 3.75, g Q and may be evaluated by the division of Equations 
(42a) and (42b) of ^ 3 ^ respectively for ri •= n— by Equation (44a) of^ 3 ^ 

cl 

— r — -v , 

for h = n. For n— > 3.75, g and g may be evaluated by the division of 
& o * 

[3] — 27 t 

Equations (44a) and (44b) of respectively for n = n— by Equation (44b) 
[3] — - 

of for n = n. In this case g Q and g^ can be expressed more simply as 


c , r, -n(l - — ) 
S (n— ) e a 

o a 


J S (n) 
a o 


(79a) 


c , r< -n(l - — ) 
S (n— ) e a 

* d 

g ‘ = „£>* S (n) 

3i O 


(79b) 


Since the computed core values to be presented herein are to be 
evaluated at the interface r = a only, the final expressions needed to 
evaluate the core displacement and stress integrals at the outer radius 
for small and large values of n become, from Equations (76) and (78). 


5 n n l ~ X " + X n® ] 


(80a) 


W n = X 12* + X U C2ot g ‘ 


(80b) 


(5 r ^ n = ^ 12 ^ ~ U - A n [(n 2 + 2a)g 


(1 + a)] 


(80c) 


(5g) n = A 12 g + X u [l - a+ 2ag] 


(80d) 


( ° x >n = * X u [2 "'i 1 


<Wn = n[ ' A >! 5 ' X H (a e ' 1>] 


(80f) 




Similarly for small and large values of n, the cylinder displacements 
and stresses given by Equations (67) are written 
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(u z ) n = ^ {-A l2 (l " hg) + A n [n 2 g + r)(2ag -1)]} (8la) 

(a xx ) n = -X l2 [l - (y + n)g] -A u {y+n-Cn 2 +2a(y+n)]g} (81b) 

(5yy) n = -X 12 Cy-(i+yn)g] -A x jd+yn-Cyn 2 + 2a(i+yn)]g} (8ic) 

For evaluation at the outer radius of the cylinder Equation (68) appro- 
priately gives the quantity n. 

Equations (59), (61), (64) and (67) apply in the intermediate ranges 
of values of n. Equations (72), (76), (78) and (81) apply in the small 
and large ranges of values of n. These equations thus provide the necessary 
quantities for the evaluation of s given by Equations (70) in the typical 
integration of Equation (69). For evaluation at the outer core radius 
Equations (65) and (80) apply as the special cases of Equations (61) and 
(64); (76) and (78). 

It should be noted that for small values of n the final bracketed 

.< x 

quantity in the integrand of Equation (69) reduces to {or}. The n— term 

1 3 

applies to the core and cylinder axial displacement and core shear stress 
quantities. The multiplication of the n of ih with U n and (u z ) R must be 
done prior to evaluation to avoid division by zero in the computations. 

For a given nondimens ional displacement or stress parameter, the 
quantity s in Equation (69) war. found to remain the same sign and converge 
to zero with increasingly larger values of n. The general nature of the 
integrand of Equation (69) is effected by the periodicity produced by the 
trigonometric functions. Greater difficulty is encountered in the numerical 
integration for the band load problem than for the limiting line load 
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problem. This is due to the presence of two periodic functions and the 

greater frequencies of oscillation associated with the consideration of f 

larger axial length values. | 

3 

Typical plots of the integrand as functions of n for the radial I 

| 

deflection for a line and two widths of band load are shown in Figure 4. f 

5 

From these plots it can be seen that the contributions to the total inte- 1 

• X 2 6 •€ 

gral occur in a continuous fashion for — = 0, — =0 but in discrete § 

3d i 

X £ 

intervals between the zeroes of the trigonometric functions for — i 0 § 

2 (S ^ 

and/or — i 0. The method of arriving at a satisfactory evaluation of j 

the integral was different for these cases. Basically, however, four 

groups were formed in both cases. These gave partial contributions towards J 

the accurate, evaluation of the integrals. These were obtained by numerical 

integration using Simpson's iile^^. In the case of — = 0, — = 0 the 

3 3 

groups in Figure 5 were chosen. These are shown by Roman numerals. For 
X 26 

— t 0 and/or — i 0, the groups were formed as the contributions of sue- 
.cessive intervals between the zeroes of the trigonometric functions. These 
are shown in Figure 6. Four groups were arbitrarily chosen, in that a 
second group may sometimes be small compared to the first and yet subse- 
quent groups may contribute to the overall integral. This is shown to 

x 26 

some extent in Figure 4 for — = 0.1 and — = 0.5. 

a a 

After these groups were formed a check was made comparing the second 
group with the first , the third with the sum of the first two and the 
fourth with the sun of the first three. Tests of 3%, 2%, and 1% were made 
to establish adequate convergence. If any group was not less than these 
specified percentages, the test was rerun with the formation of additional 
groups as necessary. 

The computer program that was set up to perform the numerical integra- 
tion is included as Table 1. This particular program was used to calculate 
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the core deflection and stress distributions in the axial direction for 

different band load widths. For convenience an equilvalence table is 

included for the algebraic and computer symbols for the different variables. 

As compiled and listed by an IBM 360 Model 50 computer in Table 1, the 

equations of the Two Function approach with interface matching used herein 

are shown. This is designated Run (1) and used to obtain results for band 

load solutions and, in the limit, for the line load solution. By a simple 

change in some of the cards of the Fortran source deck, middle surface 

matching solutions can be obtained. The Two Function approach with middle 

surface matching represents an approach that should lead to results in 

[7] 

agreement with those obtained by Yao . This solution has been designated 
Run (2). The Love Function approach with middle surface matching, which 
is considered in Appendix 2, is designated Run (3). For the line load 
solution. Run (3), results obtained should correspond to those obtained 
by Yogananda^"^. 

The program as set up for an ’accuracy’ of 1% took relatively large 
amounts of computer time, especially for the larger values of 25/a. For 
the most part, computer results were obtained on a Control Data G-20 
digital computer. On this machine for the run printed in Table 1 the 
total time for 2'5/a = .5 was close to one-half hour. As rerun on the IBM 
360 Model 50 this was cut down to six minutes. It was surprising that such 
a large reduction in time occurred in that, although the 360 is acknowledged 
to be faster, almost identical results were obtained even as far as the 
necessary number of groups to satisfy the convergence criterion. 

Numerical Results 

-5 

In all the results graphed, the material parameters are G/E = 4 x 10 
and v = y = 0.3. The geometric parameter which controls the relative thick- 
ness of the cylinder was taken as a/t * 400. Although only these particular 
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parameters are used, and results at the core and cylinder outer radi' were 
obtained, the results for a wide range of values can easily be obtained. 

Table 2 presents the core displacement and stress parameter results 
of Run (l) for the band load solution with 25/a = 0.5. These are a set of 
typical results from the computer program of Table 1. Table 3 presents 
typical results for the outer surface cylinder stress parameters. Tables 
4 and 5 present the core displacement and radial and shear stress para- 
meter results for the line load solutions of the interface matching Runs 
(2) and (3). 

Figures 7 through 12 show the distribution of core displacements and 
stresses at the outer core radius for 2<5/a = 0, 0.1, 0.25 and 0.5. For 
these same band width values the outer surface cylinder axial and circum- 
ferential stress distributions are shown in Figure 13 and 14. Results 

2 <5 

for the limiting case of line load (— = 0) are shown in all these 
Figures 7 through 14. 

No measurable difference in results for the radial displacement and 
radial stress for the two different axial displacement assumptions was 
observed in the line load solutions of the different Runs (1), (2), and (3). 
The distributions of nondimens ional axial displacement and shear stress for 
both interface and middle surface matching are shown in Figures 15 and 16. 
The cylinder axial displacement results for a line load when no core is 
present previously given in Figure 3 are also shown in the axial displace- 
ment plot of Figure 15 by the dashed curves. The significance of these 
curves is considered in greater detail in the discussion of results. 

In Figures 17 through 20 a more general picture is given of how the 
displacements and stresses at the core outer radius (at the center of the 
load) decrease with increasing band width. Figures 21 and 22 show the sane 
decrease for the cylinder stresses. In both these figures a dashed curve 
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has been used to indicate the "plane strain" solution for a finite cylinder 
and core combination where the 'length' of the cylinder is equal to the 
band load 'width'. By plane strain, it is meant that the axial displace- 
ment of both the cylinder and core is zero. In the latter calculations 
the distinction between the interface and cylinder mean surface radius was 
recognized. The derivation is given in Appendix 3. 

Discussion of Results 

Equations have been derived for the stresses and displacements within 
a long cylinder encasing a core, subject to a band of compressive stress 
uniformly distributed over a finite width 26. Results that show displace- 
ment and stress distributions versus the axial coordinate for various band 
widths and versus band width size are shown graphically in a number of 
cases for which it was assumed that the material constants G/E is 4 x 10 
Poisson's ratio v = y = 0.3 and the core radius to wall thickness ratio a/t 
is 400. 

For an increasing width of band load 26 as considered in Fig .•res 7 
through 14 the decrease in magnitudes of the displacement and stress is 
apparent. The maximum axial displacement and corresponding maximum shear, 
stress necessary tc satisfy the interface boundary condition are at similar 
locations. The location of this maximum was at an axial length slightly 
greater than the distance to the edge of the band load. In accordance with 
St. Venant's hypothesis, the axial displacement at a large distance from 
the load is independent of the particular width of load. 

The radial, circumferential and axial stresses, as well as the radial 
deflection, at the outer radius of the core, all have their maximum magni- 
tudes at the center of the band width for small band widths. For wider 
band widths, however, the peaks occur at axial locations closer to the edge 
of the band. The radial tensile stress and shear stress at the cylinder 
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to core interface are particularly important because they soecify the 
bond strength required between the core and cylinder. The maximum radial 
tensile stress values as well as the maximum shear stress values occur 
at an axial location slightly beyond the edge of the load. The maximum 
radial stress tensile values are somewhat smaller than the maximum shear 
stress values. 

In Figure 13 the axial displacements at the middle and inner surfaces 
of a cylinder with no core subjected to a line load are shown in the dashed 
curves. These curves show the difference in axial displacements at these 
two surfaces due to the presence of an axial bending moment. The correspond- 
ing cylinder-core combination line load solutions given by the solid curves 
in this figure can be seen to be a natural consequence of the particular 
matched boundary conditions. For the material and geometric properties 
considered the radial displacement and stress were not values effected v y 
the different matching (interface or middle surface) boundary conditions. 

Although a similar method of solution to that used herein was used 
for the limiting case line load solution by both Yao^^ and Yogananda^®^ 
substantially different results from those herein and from each other were 
obtained in each case. For example for the same parameters as those chosen 
herein the maximum radial stress obtained by Yao and Yogananda is about 20% 
and 37% lower respectively, than the current results. The reason for this 
difference has been noted for Yao in terms of the error in the radial stress 
expression. Yogananda's equations were computer checked with the convergence 
criteria presented herein. Rather than the results he presented, the com- 
puter results were found to lead to the values obtained in the present 
investigation (with middle surface matching). Indeed, it can be shown that 
the Two Function approach used herein by Yao and the Love Function approach 
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which Yogananda used, are algebraically equivalent. This is shown in Part 
B of Appendix 2. 

It may be of interest to examine the stresses and radial displacements 
within a section which passes through the center of the externally applied t 

pressure band in order to observe the changes that would occur as the band f 

width increases while the external force is constant. Stresses and radial 
displacement at the midsection have been calculated by the use of equations j 

derived in the present paper and are shown by solid lines in Figures 17 | 

through 22 for a band width 26/a which ranges from 0 to 2. By compariso,', 
results derived in Appendix 3 for a simplified solution of a cylinder and 
core assembly under the plane str’ain assumption that the axial displace- 
ment vanishes are shown in the same figures by the dashed lines. As expected 
the greatest difference between solutions, as measured by the vertical 
intercept between the dashed and solid lines, occurs for small band widths. 

The difference decreases considerably once the band width exceeds a length 
equal to the core radius. It is interesting to note that except for the 
axial stress component, the simplified solution provides a reasonable 
approximation to stress and axial displacement at the midsection of the 
pressure band provided the band width exceeds a core radius. 

All curves of Figures 17 through 22 show the rapid decrease in stress 
and radial displacement as the band width increases from that of a ring load. 

The rapid change empharizes the advantage of designing away from very narrow 
bands of pressure. 

In conclusion it should be noted that tensile stresses can be created 
in the bond between a case and its core even when the assembly is subject 
to a compressive band of pressure. An increase in the width of a narrow 
band of pressure significantly decreases the stresses and displacements in 
the assembly even though the external resultant force is kept constant. 


**■ 





31 . 


Acknowledgments 

The authors wish to thank the staff of the Computer Operations Section 
of the Data Processing Department at Foster Wheeler Corporation and its 
Supervisor, Mr. F. Mocenigo, for carrying out the running of the lengthy 
computer programs. Grateful appreciation is also due Mrs. Marian Cocubinsky 
and Miss Julie Bryant of the Solid Mechanics Department at Foster Wheeler 
Corporation for the rough and final typing drafts of this report. Foster 
Wheeler Corporation personnel were sponsored by the Foster Wheeler Corporation. 
New York University personnel were sponsored by National Aeronautics and 
Space Administration Grant No. NGR-33-016-067. Financial support by both 
organizations is gratefully acknowledged. 


32 . 


3IBLI0GRAPHY 


1 . Love , A . E . H . , A Treatise on the Mathematical Theory of Elasticity , 

Dover Publications, Inc., New York, 4th Edition, (1927). 

2. Harguerre, K. , "Ansatze zur Losung dar Grundgleichungen der Elastizi- 
tatstheorie" , Z. angew Math. Mech., Vol. 35, No. 6/7, (1955), pp. 242-263. 

3. Melworm, R, F., Berman, I., Shaffer, B. W., "The Displacement Function 
Approach to Cylindrical Elasticity Problems", N.Y.U. Research Rpt . S-68-4 
to NASA under Grant NGR-33-016-067 , (1968). 

4. Melworm, R. F., "A Study of the Elastic Axisymmetric Behavior of 
Cylinders with Applications" , Doctoral Thesis submitted to New York 
University, (February, 1969). 

5. Timoshenko, S. P. and Woinowsky-Kriegt:' 1 , 3., Theory of Plates and Shells , 
2nd Edition, McGraw-Hill, (1959), Cahpter 15. 

6. Donnell, L. H., "A Discussion of Thin Shell Theory", Proc. 5th 
International Cong. Appl. Mech., Cambridge, Mass., (1938), pp. 66-70. 

7. Yao, J. C., "Bending Due to Ring Loading of a Cylindrical Shell with 

an Elastic Cere", Journ. Appl. Mech., Vol. 32, (March, 1965 ) , pu. 99-103. 

8. Yogananda, C. V., "Ring Loading of a Long Thin, Circular Cylindrical 

Shell Enclosing a Soft, Solid Core: A Recalculation by the Love Function 

Method of Elasticity", Int. J. Mech. Sci., Pergamon Press Ltd., Vol. 8, 
(1966), pp. 751-757. 

9. Sneddon, I. N., Fourier Transforms , International Series of Pure and 
Applied Math., McGraw-Hill Book Co., Inc., New York, (1951), p. 11, 

Theorem 4. 

10. Kreyszig, E., Advanced Engineering Mathematic s, John Wiley and Sons, Inc., 
New York, (1962), pp. 502-508. 

11. Kaplan, W. , Advanced Calculus , Addison-Wesley , (1952), pp. 167-175. 

12. Foppl, A. and Foppl, L. , Drang und Zwang , Munchen and Berlin, 4th 
Edition, (1944). 


! 




33 . 


APPENDIX 1 


Axial Displacement of a Cylinder Without a Core 


The classical equations of thin shell theory are given by Equations 
(31). When a cylinder is subject to a line load as shown in the sketch 
below 




the boundary conditions and symmetry conditions, to be used in conjunction 
with Equations (31), may be written 


[W : -0 ; tQ ] 

x=0 x=0 


£ . limit w = 0 

2 ’ x 
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In addition the unit surface loadings p and p vanish. Equation Ola) 

X z 


can then be written 


Ru , = -yw 

x 
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(1.2a) 
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Substitution of Equation (l.2a) into Equation (31b) yields the familiar 
equation for axisymmetric cylindrical shell problems, namely 

w, + 48^ = 0 (1.3) 

’xxxx 

where the parameter 
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The solution for the radial displacement can easily be obtained from 

Equations (1.1), (1.3) (27c), and (29c). A solution is given, for example 

r 5 ] 

on pages 471-474 of Timoshenko and Woinowsky-Krieger and can be written 
in the current notation as 

w = - ggfjj e ^ X [cos Bx + sin Bx] (1.5) 

With the substitution of Equation (1.4) into Equation (1.2b), performance 
of the indicated integration, and the use of the boundary condition that 
the axial displacement at the origin of the axial coordinate vanishes, the 
expression for axial displacement may be written as 

u = [1 - e cos Bx] ( 1 . 6 ) 

The axial displacement at any distance from the middle surface can be 
found from Equations (23b), (1.4) and (1.6). When the dimensionless 
quantities that have been used herein are introduced and use is made of 
Equation (37), which expresses the equivalence of the inner and middle 
surface radii of the cylinder, the axial displacement of the cylinder may 
be written as 


2Gu 

q 

where. 
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x 10 3 = & x 1C 3 )(|) (y-e ' 6 I [y cos B 7 + (~) - sin ]£] } (1.7) 
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0 = 6 a = [3(l-n a ) (i) s ] 1/ ‘* (1.6) 
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The axial displacement ratio 2Gu z /q x 10 vs. the dimensionless axial 
coordinate x/a is plotted in Figure 3 for values of 2z/t = -1 corresponding 
to the interface (no core) and 2z/t = 0 for the middle surface (no core). It 
is of interest to note that as x/a becomes infinitely large the axial dis- 
placement ratio is given as 
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x 10 3 


(J x 10 3 ) (|) u 


(1.9) 


and is independent of 2z/t. The asymptotic value of 4.8 is reached for 
the particular parameters considered herein as shown in Figure 3. 
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APPENDIX 2 - PART A 


LOVE FUNCTION APPROACH - METHOD OF SOLUTION 


Equations (9) express the elasticity equilibrium equations in terms 
of the displacement components, W and U. The Love function X = X( r > x ) may 
be introduced in the following manner^ 


li. 

2G Dr9x 


(2A.la) 
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Substitution of Equations (2A.1) into the stress-displacement relations 
of Equations (7) then relate the radial and shear stress to the Love Function 


a = _i o ifx + (1 _ “>(i!x + 1 + if* )] 

°v. 9x L 9r^ O n ar >2 + r> + 3v2 > J 
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a 3 ; x a ,3 ; x , i 3x , 3 ; x n 
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(2A.2b) 


Substitution of Equations (2A.1) into Equations (9) leads t the 
following differential equation which governs the Love function 


32x L.22L 2 9 3 X 9 4 X 9^X 

ar 4 r 9r 3 r^" 9r 5 + r 3 9r + r 9r9x 2 + “9r 2 9x 2 + 9x 4 


(2A.3) 


The solution of this equation may be found by the use of a separation of 

[ 12 ] 

the variables technique twice. This is done m Foppl and Foppl . The 
solution may be written 


X = a 2 — , [A'l (r£) + B'n “ I,(rr-)] sin (n— ) dn (2A.4) 

A n J na a 1 a a 


where A' and B' are functions of n analogous to the functions A and B 


introduced into Equations (42) of the text; I and Ij are modified Bessel 
functions of the first kind of order zero and order one respectively. 
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Substitution of Equation (2A.4) into Equations (2A.1) and (2A.2) 
expresses the axial and radial displacements and the radial and shear 
stress components in terms of A' and B' as 
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(2A.5a) 


2GW = - — [A 1 I.(r^) + B* r£- I (i£)] cos (r£) dn 
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-r£ I.(r£)]} cos (n— ) dn (2A.5c) 
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(A 1 I , C n — ) + B'[n— I (n— ) + al (n— )]} sin (n— ) dn 
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(2A.5d’> 


Substitution of Equations (32), (34a), (35), (37), (49), (51) and (2A.5) 
each evaluated at the outer core radius, into Equations (31) leads to 
two integral equations of the form given by Equations (53). In the pre- 
sent case the A ; s have revised definitions now. Denoted by primes. 


A' 

1 1 
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(2A.6a) 


A* = ( 2a-g ) I (n) + [n 2 + a(y-g)] (2A.6b) 

1 2 w o n 

A ^ j = -g I 0 (n) + (1 + kn 4 + g -y) --jp— (2A.6c) 

A j 2 = [1 + kn 4 + g - y-a(g+y)] I Q (n) - ng I Q (n) (2A.6d) 

Expressions identical to Equations (57) through (60) and (63) hold for 

the primed quantities A', B' , A''s and for D'^ equal to the determinant 

of the A' quantities and where the terms U , W , (o ) and (t ) now 
^ n n r n rx n 
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denoted by primes have the following definitions: 

U' = - {-A' 1 (n-) + A' [2al (n-) + n- I (r£)]} (2A.7a) 

nn 120 a 1 1 oa a 1 a 


I , ( n— ) 

w ' = A ! 2 — - — § A* - I (n— ) 

n 12 n 1 1 a o a 


(2A.7b) 


I (n-) 

( VA = Kz tl 0 (nf) - -V- ] + C(l-a) ! 0 (nf)- n| I/nf)] (2A.7c) 

n— 
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(x V = -A' I .(nil) + A* [i£ I (r£) + a I,(n^)] (2A.7d) 

rxn 1 2 1 a i i aoa 1 a 


At this point, it is possible to compare the algebraic expressions 

obtained from the Two Function approach and the Love Function approach. 

This is done in Part B of this Appendix. That they are equivalent for 

the axial displacement parameter is shown for middle surface matching 

(i.e., it is shown that U ,'J, = U'/Dl ). This can be shown in a similar 

n A n A 

fashion for the other nondimens ional quantities. Rather than show this 
for the other nondimensional quantities, the numerical results, which are 
easily obtained for the different approaches have been used as a means of 
verification. Since numerical results are obtained at the core outer radius 
Equations (2A.7) reduce to 


U' = - {-A' I (n) + A' [2al (n) + n I (n )] } 
nn 120 ii o i 


(2A.8a) 


w ; = KzLdni - V n) 


(2A.8b) 
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(o )' = A’ [I (n) - — ] + A' [(1-a) I (n) -nI.(n)J 
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A similar numerical integration difficulty is encountered in dealing 
with small and large values of n in the present Love function approach as 
was considered in the Two Function approach. In this case the difficulty 
can be overcome by dividing numerator and denominator of the equivalent 
s quantity of Equation (70) by I o 2 (n). The denominator becomes 


D' = 


A' A' - A' A' 


A I z (n) 11 22 1221 
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(2A.9) 


where the functions A' , A’ , A' , and A' are defined 
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The final expressions evaluated at r = a for the numerator in the small 
and large range of n are given by 
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Equations (2A.6) and (2A.7) in the intermediate range of n and Equations 
(2A.10) and (2A.ll) in the small and large range of n provide the necessary 
quantities for the evaluation of the typical equivalent nondimensional 
quantity S in Equations (69) for the Love Fucntion approach. 
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APPENDIX 2 - PART B 

Love Function Approach - Check of Equivalence 

It is the purpose of the present section of Appendix 2 tc show, by 
examination of equations, the equivalence of the final expressions for 
the core axial displacement parameters of the Two Function and the Love 
Function approach. The middle surface matching solution is considered 
for both approaches . 

2GU 

The axial displacement parameter is given for the Two Function 

approach by Equation (60a). For the Love Function approach it may be 
written as 

00 U* 

rr S sin («”) dn (2B.1) 

o U X a 

where IP is given by Equation (2A.7a) and 



D{ = X’ X* - X' X’ (2B.2) 

X U 22 12 21 

where the X f, s are given by Equations (2A.6). 

To show the equivalence of the two approaches it is sufficient to 

show that 


U U’ 


n n 



(2B.3) 


where U and D, are given by Equations (61a) and (59) respectively. To 
n A 

do this, each side of this equation* numerator and denominator separately, 
will be put into slightly different form and it will be shown that 

U D, 

U’ = -4 ; D{ = -4 (2B.4b) 

n n 3 X n 3 


and therefore Equation (2B.3) is satisfied. 

For convenience, let the following definitions for the quantitis h } 
and h 2 be made 
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hj = n 2 + a(y + g) ; h 2 = 1 + kn“* + g-y (2B. 5a,2B. 5b) 


Then Ea”j . .ons (56) may be written 

A n = -n 2 I Q (n) + ng I^n) (2B.6a) 

A i2 = ng I o (n) " h iV n) (2B.6b) 

X = n 2 g I (n) - h n I (n) (2B.6c) 

2 1 O 2 1 

X 22 = -[h 2 -a(y-g)]n I Q (n) + [n 2 g + 2ah 2 ] I (n) (2B.6d) 


Substitution of Equations (2B.6) into Equation (59), followed by rearranging 
and reducing terms lead to 


D, = [h -a(y-g) -g 2 ] n 3 I 2 (n) + [n 2 g 2 + 2 ah g-h h ] nl %) 

A 2 O 2 12 1 

+ 2an 2 [yg - h 2 ] I (n) I (n) (2B.7) 

Substitution of Equations (2B.5) into Equations (2A.6) yields 


A' = I (n) - & I (n) (2B.8a) 

n o n i 

I. (n) 

A' = (2a-g) I (n) + [n 2 + a(y-g)] — —* (2B.8b) 

12 o n 

h 2 

A' = -gl (n) + — I (n) (2B.8c) 

21 o n i 

A 22 = [h -a(g + y) I Q (n) - ng I^(n) (2B.8d) 


Substitution of these equations into Equation (2B.2) verifies that 
Equation (2B.4b) is true. Equation (2A.7a) may be rewritten 


- I (r£) 
a i a 


From Equations (2B.6a), (2B.6b) and (2B.8a), (2B.8b) is is easily shown 


that 


- <( « i m I' 1 strewn*! 1 ■ .. - «• 
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-A* + 2aA' A nA 

U. = _ n ; A' U- (2B.10a.2B.lfo) 

n n 3 * u n J 

Therefore, from Equations (61a), (2B.9) and (2B.10) it follows that 
Equation (2B.4a) is true. Thus, the algebraic equivalence for the two 
different approaches for the axial displacement parameter is verified. 




44 . 


APPENDIX 3 


Cylinder - Core Assembly for which Axial Displacement Vanishes 


Let us analyze a core and cylinder assembly subject to a uniform pressure 
over its entire length 26 as shown in sketch (a) below. The problem is 



to be analyzed with the use of the plane strain assumption that the axial 
displacement vanishes. Hence for brevity the solution will be referred to 
as a plane strain solution, even though the axial dimension 26 need not 
necessarily be large compared to the diametral dimension 2a. The assembly 
is separated into the tv o free bodies shown in Sketch (b) and Sketch (c). 

The interface pressure between both bodies designated p^ is initially unknown. 
It will be determined from the compatibility requirement between the radial 
displacement of the core and that of the cylinder. 

T e solutions to the problems depicted in Sketches (b) and (c) can be 
taken directly from any well known text on elasticity. In terms of the 
notations used in the [ resent paper, the stresses and radial displacement 
of the core is given as 

°r = ' P I * a e = " P I 1 a x = " 2Vp I (3. la- 3.1c) 

rp I 

W=- [(1-2V)(1+V)] (3. Id) 

c 

The solution for the cylinder was found by Lame and may be found on page 144 
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[ 22 ] 

of Love . The stresses and radial displacements are found to be 

a 2 b 2 (p-p T ) L p J a 2 -pb 2 

_ o 


zz (b 2 -a 2 ) r 2 b 2 -a 2 


(3.2a) 


aVCp-pj) ± p- [ a 2 -pb 






yy (b 2 -a 2 ) r 2 b 2 ~a 2 


(3.2b) 


p x a 2 -pb 2 

°xx ~ 2lJ b 2 -a 2 


(3.2c) 


w - 


r l+u ( PrP )a2b2 P! a2 -P b2 
E [ (b 5 -a r ) '' f (1-2U)U +U )I 


(3. 2d) 


where b, is the outer radius of the cylinder. 

Equating radial displacements at the interface 


W 


r=a 


= W = w 
a 


= w 


r=a 


(3.3) 


the interface pressure p T is found to be given by the expression 



where. 




2K 2 (1-U 2 ) 

(l+u)K 2 +(l-2p)(l+y) + 


( l-2v)(K*-l) 

# 


(3.4) 


Ir. order to compare the present solution with that presented in the 
uiain section of this report, let us introduce the resultant force per 
unti circumferential q. The term q = 2fip, as given by Equation (65). 
Hence the stress components of Equations (3.1) may be written 


aa a a 

r _ _ i 

q " , 
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(3.5a-3.5c) 
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The radial deflection at the cylinder middle surface, is then given from 
Equation (3. 2d) as 


(|) 4(l+y)(~ - 1 )K 2 +(— -K 2 )(l-?y)(l+q)(KTl) 2 

= A [ 2 2 ] (3-6) 

a+b (2-) (K+1) 2 {K-1) 


The difference between this radial deflection and that at the interface 
is very small. Equations (3.5) and (3.6) were used to obtain the dashed 
curves shown in Figures 17 through 20. 

In a similar fashion the nondimensional axial and circumferential 
stress ratios can be expressed in terms of nondimensional quantities 
introduced herein. From Equations (3.2b) and (3.2c) one may write 


aa 


xx 


— - K 

JL- [ 2 

A K “ 1 

a 
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-k 2 [i - P] 

u7% 

[i + L. ] 2 

2(f) 


K 2 ) 


(2-) (K 2 -l ) 
a 


(3.7a) 


(3.7b) 


These equations were used to find the dashed curves shown in Figures 22 
and 23. 
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TABLE 1 

COMPUTER PROGRAM TO CALCULATE AXIAL VARIATION OF CORE DEFLECTION AND 
STRESS PARAMETERS FOR DIFFERENT BAND LOAD WIDTHS 


The program is given on the following pages. The algebraic and 
computer nomenclature for the significant variables in the problem 


considered are: 

Algebraic 

TT 

G/E 

V 

V 

a/t 
25/a 
6 /a 
a 

g 

k 

x/a 

n 

I (n) 
o 

Ij(n)/n 

I x (n) 

S (n) 
SjCn) 
g 


Computer 

PI 

GE 

VF 

V 

AT 

DL2 

DL 

AL 

G 

PK 

XX 
X 

XIO 

XI1N 

XII 

XXO 

XXI 
GS 
FI 
F2 
F3 
F4 


TABLE 1 (continued) 48 _ 

DISK OPERATING SYSTEM/360 FORTRAN 360N-FO-451 21 

C REQUEST FOR STORAGE 

DIMENSION U( 6) ,M(6) ,FF1 (44) , FT 2 ( 44) t FT 3(44) ,FF4(44) , FF5( A A ) , 

1 ‘ FF6( 44) ,SP(6) ,XL (6) ,XD(6) ,S(4,6i ,A(4,6),8(4,6) 

C CONSTANTS 

PI=3.141592654C0 

C INPUT MATERIAL AND SIJE PARAMETERS 
10 REA!) ( 1,20) GE v V , VF,AT,DL2 
20 FORMAT (5F10.0) 

C CALCULATIONS OF OTHER MATERIAL AND SIZE PARAMETERS 
AL=2.*t 1-VF ) 

G=V-2.*( l.-V*V)*GE*AT 
PK= 1 . /( 12«*AT**2« ) 

DL = DL2/2. 

C RUN NUMBER 
N = 1 
C HEADING 

WRITE (3,30)N,DL2,GE,AT,VF,V 

30 FORMAT (1HL , 57X, • TABLE 2 • / / 37X , • E L AST I C ANALYSE OF CYLINDER WITH 
1COP.C - RUN ( »I 1 »•) '/MIX. 'UNIFORM BAND LOAD OF * . DTH 20/ A= • F 5 • 2 / / ' 
1 CX » * G/E ='F12.8,3X,*A/T ='F7.2»3X»' POISSON RATIO - CCRE =»F6.3/77/ 
1 , 'CYLINDER =*F6.3/ ' 1 SX , » X/ A • 2X , » 2 GU/OUO 3 • 2X , • M • 5X , • 2CW/C • 3X , • M ' * 
IX, »TRXA/0'3X. 'M'5X , 'SRA/Q*3X, 'K'5X* ' SXA/C 3X f 'M' 5X. * STA/0' 3X, 'M 1 / ) 
C OUTPUT LOCATIONS AS DEPENDENT UPON D/A 
IF (DL) AO, 50,40 
40 IT = 19 
GO TO 60 
5C IT = 21 
60 DO 1230 I = 1 , I T 
IF (1-1) 70,70,30 
70 KS = 0 
30 XI = I 

IF (DL) 100,90,100 
90 XX = ( X I -1 - 1/100. 

GO TO 160 

100 IF (l-l) 130,130,110 
110 IF (1-5) 140,140,120 
120 IF (1-15) 150,150,140 
130 XX * 0 

GO TO 160 
140 XX = XX+DL/5. 

GO TU 160 
150 XX * XX4-QL/25. 

C SIX VALUES FOR DEFLECTIONS AND STRESSES DETERMINED 
160 DO 1260 L.S = 1,6 

IF ILS-2) 200,1 ZG » 1 30 
170 L a 3 

GO TO 210 

130 IF (LS-4) 190,200,200 
190 L « 2 

GO TO 210 
200 L = LS 
210 M. ( L ) - 0 


TABLE 1 (continued) 
10/09/67 FOKTMAIN 

IF (KS1230, 220,230 
22 0 SP l L 1 - n 
230 KK=C 

A(1«L)-C 
Kl=2 
JI = 1 

C SLLECTION flF NUMERICAL INTEGRATION LIMITS 
240 00 250 Jl = J I «4 
250 S(Jl,L)=C 

IF (XX) 450,260,450 
260 IF (L-l 1270, 1260,270 
270 IF (L-3)280, 1260,280 
280 IF ( DL 1 370 » 290 , 370 
C NO PERIODIC I TY -DL=0,XX=0 
290 IF (KK 1310,300,310 
300 PI 1,L)=50. 

310 DO 360 K1=KI,4 

IF (KK) 320,340,320 7 

320 IF (Kl— Kl 1340,330,340 
330 KT=4 

GO TO 350 
340 KT=Kl-l 
350 A(K1,L)=B(KT,L) 

360 B( K J. ,L ) = A ( Kl , L ) ♦5C« : 

ko- l ,-v' 

GO TO 810 ; 7 Jl 7_ _ 

C ZEROES OF SIN (X*DL) - XX=0 
370 IF (KK 1390,380,390 

380 B( l ,L 1 -PI /OL 

390 DO 440 K3=K I ,4 ' . : 

IF (KK 140 0,420, 400 ' - •> 

40C IF (K3-KI 1420,410,420 
410 K T =4 

GO TO 430 
420 KT -K3— l 
430 A( K 3 , 1. 1 =8 l KT ,L) 

440 P(K3,L)=A(K3,L)*PI/DL 

KO= 1 _ J 

GO TO 810 

C ZEROES OF COS (X*XX) OR SIN (X*XX) - DL=C 
450 IF (OL 1600,460,600 
460 IF (KK1510, 470,510 
470 IF (L-l) 4 30,490,480 
480 IF (L - 3 150 0,490, SCO . .. 

490 B(1,L)=PI/XX 
GO TO 510 

500 BC1,L)=PI/C2- > 

510 DO 560 K 2 = K I ,4 

IF (KK 1520,540,520 
520 IF (K2-KI 1 540 ,5 30 , 5*t0 _ 

530 KT=4 

GO TO 550 

540 KT=K2-l ... 


(continued) 
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TAP” 

10/09/67 FORTMAIN 

550 A(K2,L)=B(KT,l ) 

560 BlK2,L)=A(K2,L)*PI/XX 
IF (L-l 1570,530,570 
570 IF (L- 3)5 90 ,580,590 
580 K0=? 

Gfj TO 810 


590 KO=l 

GO TO 8 1C 

C ZEROES OF BOTH COS (X*XX) OR SIM (X*XX) AMD SIN 
60 v IF ( KK ) 700,6 10,700 
610 Xl(L)=Pl/XX 
XD { L ) =?I /DL 
IF (L-l )620,640,620 
620 IF (1-3)630,640,630 
630 XL ( L )-XL( L ) /2. 

640 IF ( XL(L)-XD(L) )65C,670,680 

650 IF (A3S(XL(L)-XD(L) )-.Cl ) 670,670,660 

660 PI 1 ,L) = XL(L) 

XL(L)=XL(L)+PI/XX 
GO TO 700 
670 XD(L)=2.*XD(L) 

GO TO 660 

680 IF ( ABS(XL(L)-XD(L ) )-.01 ) 670,670,o90 
690 R( l ,1 ) =XD( L ) 

XD(L)=?.*XD(L) 

700 DO 809 K4=K 1,4 

IF (KK) 710,730,710 
710 IF ( K4— K I )730, 720,730 
720 KT=4 

GO TO 740 
730 KT=K4-i 
740 A(K4,L)=B(KT,L) 

IF ( XL(L )-XD(L ) )750, 770,780 
750 IF ( A8S(Xl(L)-XD(l.) ) - .01 ) 770,770,760 
760 R ( K 4 , L ) = XL ( L ) 

XL(L)-XL(L)+PI/XX 
GO TO 3C0 

770 Xi> ( L ) = XD( L ) + P I /DL . 

GO TO 760 ** 

780 IF (ABS (XL(L)-XD(L) ) - • 0 l ) 770,770,790 
790 R(K4,l )= XD(L) 

XD(L)=XD(L)+PI/DL 
800 CfJNTI MUE 

KO = 2 . 

810 Ml =0 

C FORMATION OF FOUR GROUPS FOR CONVERGENCE TEST 
DO 1170 J4 = JI ,4 
M ( L ) = M f L ) ♦ l 

IF ( M ( L ) - 98 ) 820,820,1260 
820 no 1170 J * 1 v 11 
M 1 = M 1 ♦ 1 
IGN=(-1 ) **J 


IF ( J- 1 ) 830,830,840 


(X*DL ) 


TABLE 1 (continued) 

10/09/6 7 FORTMAIN 

650 A ( K? » L ) =3 ( K T , L ) 

560 b(K2,L)=A(K2,L)+PI/XX 
IF (L-l >570,580,570 
570 IF (1-1)590,580,590 
580 KO = ? 

Gf J TO RIO 
590 K0=1 

GO TO 810 

C ZEROES OF BOTH COS IX*XX) OR SIN (X*XX) AND 
600 IF (KK) 700,610,700 
610 XL ( L ) = p I /XX 
XOIL)= p I/DL 
IF (L-l >620,640,620 
620 IF IL-3)630,640,630 
630 XL( L)-XL(L)/2. 

640 IF (XL(L)-XU(L) >650,670,680 

650 IF ( A3S(XL(L»-XD(L) )-.Cl ) 670,670,660 

660 P ( l,l)=XLU) 

XL(L)=XL(L>+Pl/XX 
GO TO 700 
670 XDIL)=2.*XD(L> 

GO TO 660 

680 IF I ABS(XL(L)-XD(L) ) — - 0 1 ) 670,670,690 
690 PC I ,L) =XD(L> 

XO<L)=?.*XD(L) 

700 no 800 K4=K I ,4 

IF (KK) 710,730,710 
710 IF (K4-KI >730,720,730 
720 KT=4 

GO TO 74C 
730 KT=K4- 1 
740 A(K4,L)-B(KT,L) 

IF ( XL(L)-XDIL) >750, 770,780 
750 IF (ASS(XL(L)-XD(L) ) — • C 1 ) 770,770,760 
760 R I K4 , L ) = XL I L ) 

XLIL)=XL(L)+PI/XX 
GO TO SCO 

770 XD(L)=Xn(l.) + PI/DL 
GO TO 760 

780 IF I A3S (XL(L)-XD(L) ) - . 0 1 > 770,77C,790 
790 P.(K4,L) = XO(L) 

XO(L)=XD(L)+PI/OL 
800 CONTINUE 
K0 = ? 

810 M 1=0 

C FORMATION OF FOUR GROUPS FOR CONVERGENCE TEST 
DO 1170 J4 = JI ,4 
M(L)=M(L)+l 

IF ( MI L ) -98 ) 820,820,1260 

820 DO 1170 J * 1 ,11 . 

M1 = MIM 
I GN = I — 1. ) **J 

.. IF ( J-l ) 830,830,840 


TABLE 1 (continued) 5 o 

10/09/67 FORTMAIN 

830 X*A(.inU 
GG - 1 . 

GO TO 890 

840 IF (J-ll) 860,850,850 
. 850 X-B ( J4 » L ) 

GG=l. 

GO TO 890 

860 X=X+IB( J4,L)-A( J4,L) ) /10 . 

IF (IGM) 870,870,880 
870 GG = 2. 

GO TO 390 
880 GG = 4. 

C SETUP FOP STURAGE 

890 IF (KK) 930,900,930 
900 IF IKU-l) 920,910,920 
' 910 IF ( L-2 ) 930,930,1 130 
920 IF ( L— l ) 930,930,1130 
C LOAD ' 

930 IF (DL) 940,940,950 
940 F5={G-V)/PI 
GO TO 970 

950 IF (X) 960,940,960 

960 F5= (G-V ) /PI *SIN(X*DL)/(X*DL) 

C MODIFIED" CESSEL FUNCTIONS " ' 

970 T = X/3. 75 

IF (T) 990,980,990 

980 <10 = 1. ' 

XI IN = .5 
XU = 0 
GO TU 1040 

99-0 IF C T- 1 . ) 1000,1000,1010 

1000 XI 0= l. * 3. 51 56 2?9E0*T*T+3.C899424EC*T* *4+1.20 67492 E0*T**6+. 2659 732 
1 T**8+. 9360766 * T** 1 C+. 004 581 3* T** 12 
XI 1N= { l./2.+.87o9C594EC*T*T+.5149S369* t ** 4 + . 1 5C 3 49 34EC * I **6 + 

1.0265873 3 EO*T**3 + .O33O1532E0*T **10 + . CO'- >24 1 1 EO* T** 1 2 ) 

XIl*XllN*X 
GO TO 1030 

...1010 XX0 = I. 39394228 E0+. 01 328 592 EG*T*=* (>-l )+.00225319Eu*T**{-2)- 

i.00157565E0*T**(-3)+ .009 1628 1EC*T ** f -4 }-. 020 57706E 0* T** ( -5 ) + 
1.02635537E0*T**(-6 01647633 E0*T**( -7 )+.00392377t0*T**(-6 ) ) 

XX 1 = C .39894228EG-.03988024E0*T**(-l )-. 003620 1 8£0*T** { -2 ) + 

1 .OOl638ClEO*T**(- 3 )- .0 1C 31 555EQ*T ** ( -4 ) + .02282967c 0*T** { -5 ) ~ 
l.02695312E0*T*< ‘ -6 ) ♦ . 0 1 787654E0*! ** ( - 7 )- .004 2CC 59cO*T*< (-8) ) 

IF IX-50.) 1020,1020,1050 
1020 XI0-XX9*EXP(X)/S0RT(X) 

XI l=XXl*cXP(X)/SGRT(X) 

XI1N=XI1/X 

1030 IF I X— - 001 ) 1040,1040,1090 
C SMALL AND LARGE VALUES OF X 
1040 GS=XI IN/XIO 
GO TO 1 C60 

1050 GS = XX1/(X*XX0) 

1060 FI = -l.+(G+X*X/(2.*AT) ) * 3 S 
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TABLE 1 (continued) 

10/09/67 FORT M A I N 

F? = G+ X*X/ ( 2 . * A T > - ( X* X* ( 1 . ♦ AL / AT) ♦ AL * < V + G ) ) * GS 
F 3 = G-{l# + V*X*X/(2.< t AT) + p K< : X**4 + G-V)vOS 
FA = -.( l. + V*X*X/(2.<‘AT) + PK*X**4- ( 1 . ♦ AL ) < { N/-G ) ) ♦ ( { X* X + ? . */.L X=C + ? . 1 

? *Al*( l.-V*PK*X**4)+V*X*X*AL/AT )*GS 1 

D=F ! *F 4-F2 ^F 3 1 

Sl = (-F2*(GS-i. )-Fl*( (X*X+2.*AL)*GS-( 1 . ♦ AL ) ) )*F 5/D j 

Wl=(P2*SS*Fl*{2.SAL<GS-i. > )*F5/D I 

SX1=(-F2*F l*( 2.-AL+X*X*GS ))*F5/D | 

ST1*(F?*GS+F1*( l.-AL+2.*AL*GS) )*F5/U j 

IF (X-25.) 1070*1080. 10B0 ! 

1 0 7 C Ul = (-r < <-Fl*X*X*CS) *F5/D I 

T 1 = X*X*{-F2*GS-F1*< AL*GS-1 . ) >*F5/D 

GO TO HOC j 

1080 Ul = (-F2 + Fl*X 1 < t X£GS)*F5 / ( X * D ) ] 

T1 = X*(-F2*GS-F1*< AL*G5-1. ) I*F5/D | 

GO TO 1100 1 

C IN TERM EDI ATE VALUES OF X I 

1090 Fl=-X*X*XIO+( X*C«-X**3/(2. *AT> J * X I 1 I 

F2-( X*G + X**3/ (2.*AT) ) *X I 0- ( X *X* ( l . + AL / AT ) + AL* ( V + G ) ) * X I 1 J 

F 3 = X^X^G^XIO-XM l.+V&X+X/l 2 • * AT ) fPK*X**A+G-V)*XI 1 J 

fa = -X*( l.+V*X*X/(2.<AT)+PK*X**4-( l.+AL)*(V-G) )*X.IC + ( f X< : X + 2. * AL ) 
1*G+2.*AL*( l.-V+PK*X**4)+V*X*X*AL/AT )<=XI 1 
D=F1*F4-F2*F3 

Sl=(-F?*IX*XI l-X*X$XIC)-Fl*( (X*X + 2.*AL)*XI l-( 1 . + AL I * X* X I 0 ) ) *F 5/D 
Ul=X< ! (-F2 <£ XI0 + Fl< ! XI 1 ) * F 5 / D 
W1={F2*X$XI 1+F1=M2»*AL*XI 1 -X*X 10 ) )*F5/D 
Tl=(-F2*X*X*XI 1-F1*( AL*X*X I i-X*X*X!G ) )*F5/D 
SXl={-F2*X*X*XI0+Fl*n2.-AL)*X*XI0+X*X*Xll))*F5/0 
STl = IF2< i XIl^X+Fl<=( U.-AL)*X*XlO + 2.*AL*XI 1 ))*F5/U 
C STORAGE 

1100 IF C Kk. ) 1110,1110,1120 
1110 F F 1 ( H 1 ) = U 1 
F F 2 ( M l ) = W 1 
F F 3 ( M 1 ) = T 1 
FF4IM1 ) =S 1 
FF5IMI ) = SX1 
FF6IM1 )=STl 

1120 IF (X-. C01) 1140,1140,1150 
1130 Ul=FFHVl) 

W1=FF2(M1) 

Tl =FF 3 ( . v 1 ) I 

Sl=FF4(yi) j 

SX 1 = FF 5 ( Ml ) I 

ST l = F F 6 ( Ml ) I 

IF (X-. 001) 1140,1140,1150 j 

1140 UN=U1*XX J 

WN = ln l I 

TM^TnXX 1 

SN -S 1 ! 

SXN = SX l . ... . f 

STN= S T 1 
GO TO 1160 

1150 UN=Ul*SIN! X*XX) ... ... I 
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TABLE 1 (continued) 

10/09/67 FORTMAIN 

WN = W1*C0S< X*XX ) 

TN=T1*SIN(X*XX) 

SN=Sl*COS(X*XX) 

SXN*SK1*C0S( X*XX) 

STN=STl*COS( X*XX) 

1160 U< 1)=UN*ICOO. 

UI2)=WN 
U(3)=TM 
U{ 4)=SN 
U( 5>=SXN 
U( 6)=STN 

1 1 70 S< J4,L)=SM4 f L)+GG*U(L)*(P( J4.M-M J4,L))/3G. 

C CONVERGENCE TEST 

DO 1180 II * 1,4 

IF (ABS(S(II»l))-*010*ABS!SPlL)))1180tll80,1190 
1180 CONTINUE 

S(1,L)=S(1,L)+S;2,L)+S(3,L)+S(4,L) 

CO TO 1260 

1190 DO 1210 LI = 1,3 
XL 1 =1 1 

IF (AQSISILI+ l,L))-I .04-.01*XLl)*AB$(S(l,L) )) 1200 , 1200, 1220 
1200 S( 1,L) = S(1 ,L) + S(L1U*L) 

1210 CONTINUE 

GO TO 1260 ' ' 7 

1220 S( l ,L ) -S ( l ,L ) +S( l 1 +1 , L ) 

KK = 1 

IF (L 1-2) 1230, 1240, 1250 

1230 S(2,L)=S(3,L> 

S(3tL) s S(4,L) 

Kl *4 
JI = 4 

GO TO 240 
1240 S(2,L)=S(4,L) 

KI =3 
J l = 3 

GO TO 240 
1250 K I =2 
Jl=? 

GO TO 240 

C END CONVERGENCE TEST 
1260 CONTINUE 

DO 12/0 I K = \, 6 

127C SP( IK)=AMAX1(ABS(S(1,I<) >,SP( IK) ) 

K$=l 

C PRINTOUT OF RESULTS 

1280 WRITE ( 3,1290)XX,( (SI 1 ,N1 ) ,M(Nl) ) » Ni = 1 , 6 ) 

1290 FORMAT ( 1 5X , F 7. 3 , 6 I E 1 1 - 3 , 1 3 ) ) 

GO TO 10 
END 
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tabu 3 

CLASTIC ANALYSIS OF CYLINDER WITH CORE - RUN (I) 

UNIFORN BAND LOAD OF WIDTH 2D/A = 0.50 

0.00004000 \/T = 400.00 POISSON RATIO - CORF = 0.300 

CYLINDER = 0.300 

CYLINDER STRESSES AT l BAR = I.OC 


AXIAL CIRCUMFERENTIAL 


x/\ 

SX K* \/G 

M 

SYY*A/C: 

9 

0.0 

-C.333E 

01 

51 

-0. 745E 

0 3 

7 

0.050 

-C.303E 

01 

40 

-C.740F. 

03 

n 

C.1C0 

-C.299C 

Cl 

51 

-C.733F 

C 3 

10 

c. tsc 

-C.241C 

02 

36 

-0. 729E 

0 3 

9 

0.200 

-c • i ?y£ 

03 

22 

-C.7C8E 

03 

10 

C. 250 

C. 129E 

01 

22 

-0.378E 

C 3 

13 

0. 3 CO 

0. 129C 

03 

22 

-0.560E 

C2 

22 

0.350 

C.254E 

02 

46 

-0.214E 

02 

34 

C .400 

0.45BE 

01 

51 

-0.169E 

02 

3 9 

•' 0.450 

0.363E 

01 

51 

-0.1C8E 

C2 

50 

0.500 

- 0.359E 

01 

51 

-0.817E 

01 

17 

0.55C 

C.312E 

01 

51 

-0. 112E 

Cl 

ie 

0.600 

0.264E 

01 

51 

0.736E 

■?0 

16 

0.650 

0.3146 

01 

51 

0.A35E 

01 

16 

0.7C0 

C.298E 

01 

51 

-0.367E 

Cl 

17 

‘ 0.750' 

C.223E 

01 

51 

0.482E 

Cl 

13 

. C.80C 

C • 2 7 9 E 

01 

51 

0.421E 

Cl 

In 

0.B50 

0.2C4E 

01 

51 

-0. 121E 

Cl 

19 

6.000 

0.259C 

01 

51 

0.467F 

CL 

20 

0. 950 

C.277L 

Cl 

AO 

0. 704E 

C l 

16 

i . 000 

0.2060 

01 

46 

0 . A97E 

CO 

21 
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Figure 2. STRESS. LOAD AND BODY FORCE 

SULTANTS ACTING ON MIDDLE 
RFACE ELEMENT AND THEIR 


[•I 


INSIDE SURFACE DISPLACEMENT 


Middle surface 
DISPLACEMENT 


6/E 

= 4x10" 5 
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A 

« .3 


Figure 3. AXIAL VARIATION OF THE AXIAL 


SUBJECTED TO LINE LOAD AT 


THE INSIDE AND MIDDLE SURFACES 











Figure 4 RADIAL DEFLECTIONS AS TYPICAI 
TO BE INTEGRATED NUMERICALLY 






ZEROES OF cos (o-^-) 



TRIGONOMETRIC FUNCTIONS ARE PRESENT 
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c* 

Cu 


distributions o f radial displacements 

AT OUTER CORE RADIUS FOR ?l =Q |. 25. 5 
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Figure 9. DISTRIBUTIONS OF SHEAR STRESSES AT 
OUTER CORE RADIUS FOR 2 § = 0. I. ,25. .5 
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Figure 12. DISTRIBUTIONS OF AXIAL 

STRESSES AT OUTER CORE RADIUS 
FOR Zn s 0. .I..25..5 
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Figure |3 
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SHELL THEORY (NO CORE) 

CYLINDER (SHELL THEORY) - CORE (ELASTICITY) COMBINATION 


26 (AXIAL DISPLACEMENT) x 10* 



Figure 15. AXIAL DISPLACEMENTS FOR LINE LOAD- 
DIFFERENT MATCHING ASSUMPTIONS 
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Figure 16. SHEAR STRESS DISTRIBUTIONS 
FOR LINE LOAD 


IFFERENT MATCHING ASSUMPTIONS 
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CORE (ELASTICITY) - CYLINDER (SHELL THEORY) COMBINATION 

PLANE STRAIN SOLUTION 
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Figure 18 . CH ANGES IN RADIAL NORMAL STRESSES OF 
CORE AT BAND LOAD CENTER vs BAND LOAD 

WI DT H 
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CORE (ELASTICITY) - CYLINDER (SHELL THEORY) COMBINATION 

PLAIN STRAIN SOLUTION 



Figure 1 9. CHANGES IN CIRCUMFERENTIAL NORMAL STRESSES f 

RE AT BAND LOAD CENTER vs BAND 



CORE (ELASTICITY)- CYLINDER (SHELL THEORY) COMBINATION 
PLANE STRAIN SOLUTION 
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Figure 20. CHANGES IN AXIAL NORMAL STRESSES OF 
CORE AT BAND LOAD CENTER vs BAND 

LOAD WIDTH 
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CORE (ELASTICITY ) - CYLINDER (SHELL THEORY) COMBINATION 
PLANE STRAIN SOLUTION 
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Figure 21. CHANGES IN AXIAL CYLINDER STRESSE S 
AT BAND LOAD CENTER vs BAND LOAD 

mm 
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Figure 22. CHANGES IN CIRCUMFERENTIAL CYLINDER 


BAND LOAD WIDTH 




